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We investigate the cosmological role of a Tracking Field (f> in Extended Quintessence scenarios, 
where the dynamical vacuum energy driving the acceleration of the universe today possesses an explicit 
coupling with the Ricci scalar, R, of the form F((j>)R/2, where F(cj>) mimics General Relativity today, 
F(4>o) = 1/8ttG. We analyse explicit non-minimally coupled (NMC) models where F(cj>) — l/8nG + 
■ £(4> 2 — 4>o), with £ is the coupling constant and cj>o is the Q value today. Tracker solutions for these 

NMC models, with inverse power-law potentials, possess an initial enhancement of the scalar field 
dynamics, named i?-boost, caused by the effective potential generated by the Ricci scalar in the 
Klein-Gordon equation. During this phase the field performs a "gravitational" slow rolling until 
the true potential becomes important. We give accurate analytic formulas describing the i?-boost, 
showing that the Quintessence energy in this phase scales with the redshift z as (1 + z) 2 . When the 
' J?-boost ends, the field trajectory matches the tracker solution in minimally coupled theories. We 

compute perturbations in these Tracking Extended Quintessence models, by integrating the full set 
of equations for the evolution of linear fluctuations in scalar-tensor theories of gravity, and assuming 
Gaussian scale-invariant initial perturbations. The Integrated Sachs Wolfe (ISW) effect on the Cosmic 
Microwave Background (CMB) angular spectrum causes a change SCt/Ce — 2[1 — 8nGF((j>d e c)] at 
£ < 10, where dec stands for decoupling. Similarly, the CMB acoustic peaks multipoles shift compared 
to ordinary tracking Quintessence models by roughly an amount 81/1 ~ [8nGF(<f)dec) — l]/8. The 
turnover wavenumber ktum in the matter power spectrum shifts by an amount Skturn/kturn — [1 — 
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' 8ivGF((f) eq )]/2, where eq stands for matter-radiation equivalence. All these corrections may assume 

positive as well as negative values, depending on the sign of the NMC parameter £. We show that 
the above effects can be as large as 10 — 30% with respect to equivalent cosmological constant and 
ordinary tracking Quintessence models, respecting all the existing experimental constraints on scalar- 
tensor theories of gravity. These results demonstrate that the playground where the next decade data 
will have their impact includes the nature of the dark energy in the Universe, as well as the structure 
?—j I of the theory of gravity. 
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I. INTRODUCTION 



Following the experimental evidence in favor of an accelerating universe [Jl| — Kfj , the cosmological role of a minimally- 
coupled scalar field, considered as a "Quintessence" component (which we call here Q or <f>, equivalently) , providing 
the dominant contribution to the energy density of the Universe today in the form of dynamical vacuum energy, has 
been widely analyzed in the recent past Q - pp| . The main feature of this dynamical vacuum energy component, 
that could also allow to distinguish it from a cosmological constant, is its time-dependence and the wide range of 
possibilities for its equation of state. Moreover, by the principle of general covariance this time varying scalar field 
should also develop spatial perturbations. 

In a number of recent papers the available range of couplings between Q and other physical components has been 
explored. The possible coupling between a Quintessence field and light matter fields has been investigated in jlj] and 
is subject to restrictions coming from the observational constraints on the time variation of the constants of nature. 
Also, some recent work explores the cosmological consequences of a coupling between Quintessence and matter fields 
[Of. The possible coupling between Q and the Ricci scalar R was the subject of our recent analysis pL where the 
observable impact of such a coupling has been throughly investigated. Within this 'Extended Quintessence' (EQ) 
model two sets of theories were studied: Induced Gravity (IG, originally proposed in |l5|]) and Non-Minimal Coupling 
(NMC, sec JXg) for an extensive overview). In both models the term describing the coupling between the Ricci scalar 
R and the Quintessence field </> in the Lagrangian has the form F(tfi)R/2, where F(<fi) — £</> 2 ; the difference between 
the two models is that in the NMC model the Lagrangian contains also the ordinary gravitational term R/l6irG, that 
is instead absent in IG. Both models belong to the general class of scalar-tensor theories of gravity. The effects on 
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Cosmic Microwave Background (CMB) anisotropics, as well as on Large Scale Structure (LSS), have been described 
and numerically computed for several values of £ > 0; they turned out to be similar for IG and NMC, although 
for different values of the coupling constant. The main results of our previous study were a modification of the low 
redshift dynamics, enhancing the Integrated Sachs Wolfe effect (ISW), as well as a shift of the CMB acoustic peaks 
and matter power spectrum turnover because of the dynamics of the Hubble length. Predicting in detail the shape 
and the position of the acoustic peaks in a given cosmological model is of course important in view of the formidable 
observations of the next decade Jl7|], but it is also of great importance at present, because of the strong evidence in 
favor of the existence of subdegree acoustic oscillations in the latest CMB data |l8|,[l9| . 

As detailed in our previous paper p^j , we named our models 'Extended Quintessence , in analogy with the 'extended 
inflation' models |20| , proposed in the late eighties to rescue the original idea of inflation based on a first-order phase 
transition. In those models a Jordan-Brans-Dicke (JBD) scalar field M] was added to the action to solve the 'graceful 
exit' problem of 'old inflation'. The idea of using a non-minimally coupled scalar field as a decaying cosmological 
constant dates back to a paper by Dolgov [g2[ in 1983, which however radically differs from the present class of models 
by the absence of a true potential energy for the scalar field. The background dynamical properties of a Quintessence 
field coupled to the Ricci curvature have also been recently studied in p3j , [EJJ and |2q], |2(|] where, in particular, 
the existence and stability of cosmological scaling solutions was analyzedrModels in scalar-tensor theories of gravity 
involving scalar matter couplings have also been studied [p7| . 

Let us come now to the subject of the present work. The existence of a considerable amount of vacuum energy in the 
Universe, as observations seem to imply, brings together two conceptual problems. The first regards the magnitude 
of this energy, namely the fact that we are right now in the phase in which the vacuum energy is starting to dominate 
over matter. The second, which we name "fine-tuning" problem, which arises merely from the fact that if the vacuum 
energy is constant, like in the 'standard' cosmological constant scenario, then, at the beginning of the radiation era 
the Q energy density should have been vanishingly small compared with both radiation and matter. 

Important recent work demonstrated that the Quintessence scenario can avoid such fine tunings, while a cosmological 
constant fitting the data would be affected by both problems. Indeed, in the Quintessence scenario one can select a 
subclass of models which admit " tracking solutions" || : following early work by Ratra and Peebles |7) , it was shown 
how the observed amount of scalar field energy density today can be reached starting from a very wide set of initial 
conditions, covering tens of orders of magnitude. In particular, the Quintessence could have been initially at the level 
of the ordinary matter, thus being likely one of the products of the decay of the vacuum energy responsible for the 
inflation era. 

Here we face the extension of tracking phenomenology in EQ; we find important results, both for background 
dynamics and perturbations, and we name this scenario Tracking Extended Quintessence (TEQ). We consider a 
NMC model for the coupling between Q and R. We first analyse the background evolution, accurately studying the 
Quintessence dynamics in the radiation dominated (RDE) and matter dominated (MDE) eras. Then, we compute 
the evolution of cosmological perturbations pointing out important observable effects on CMB and LSS. 



The work is organized as follows: in Sec. II we recall the relevant equations, and we refer to previous works [JL4 10 
for a complete presentation of the dynamical system for the background and for the perturbations in non-mimmally 
coupled scalar field cosmologies; Sec. Ill contains the analysis of tracking solutions; Sees. IV and V contain the 
description and a discussion of the effects on CMB and LSS, respectively; Sec. VI discusses how these effects change 
if we vary the form of the Q-field potential. Finally, Sec. VII contains a summary of our results and some concluding 
remarks. 



II. COSMOLOGY IN SCALAR-TENSOR THEORIES OF GRAVITY 



In this Section we review the relevant equations in scalar-tensor theories of gravity, for the background evolution 
and linear perturbations. For more details, we refer to previous works ||@- 
Following , scalar- tensor theories of gravity are represented by the action 

i/(0,i?)-^(0)^0 ;M -V(0)+L /iMd , (1) 

where R is the Ricci scalar, and Lf\ U id includes all components except for 0, that we suppose coupled only with 
gravity. 

The models for f(4>, R) that we investigate here are simply a product of a function of </> and R; also we assume a 
standard form for the kinetic part: 

f{<t>,R)=F{ct>)R , w{cj>) = l, (2) 
As in ||, we assume the cosmological metric has the flat Friedmann-Robertson- Walker (FRW) form: 

ds 2 = a 2 [-dr 2 + (dij + h^dx'dx 1 ] , (3) 
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where a(r) is the scale factor, r is the conformal time and h%j represents the perturbations in the synchronous gauge, 
to be discussed later. Our units are c = 1, and we adopt the signature (— , +, +, +). 

In the next subsections we describe the evolution equations for the relevant quantities, for the background and 
linear perturbations. 



The background FRW equations read 



A. Background 



■H- = [ ~ ] = ^ Wpfimd + + a 2 V - 3HF 



(4) 



H = H 2 - ^ (a? (p fluid +Pfluid) + 4> 2 + F - 2HF 



(5) 



where the dot denotes differentiation with respect to the conformal time r; note that in our framework 1/F today 
must be 8irG. 

The Klein-Gordon equation reads 



2H(j)=—F 4> R-a^ 



(6) 



where the subscript 4> denotes differentiation w.r.t. <p. Furthermore, the continuity equations for the individual fluid 
components are not directly affected by the changes in the gravitational field equation, and for the i— th component 
we have 



pi = ~3H(pi +Pi) . 

Let us give now useful expressions for the Ricci scalar in scalar-tensor theories: 
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4> 2 [ F TiF 

-p fluid + fluid + — - 4V + 3 I + 2— 



This can also be written as a function of the Hubble parameter and its derivatives: 

R= %{H 2 + H 2 ) . 



(7) 



(8) 



(9) 



As discussed in more detail in juj , let us recall the two experimental constraints that these models have to satisfy. 
One comes from the time variation of the gravitational constant in scalar-tensor gravity theories |^| 



Gt 
G 



Ft 
F 



(10) 



where the subscript t denotes differentiation w.r.t. the cosmic time t; Gt/G is bounded by local laboratory and solar 
system experiments [E9J to be 



9* 

G 



< 10 11 per year 



(11) 



Another independent experimental constraint comes from the effects induced on photon trajectories 



1 We consider l/87rF as effective time dependent gravitational constant. We checked that, for the small values of £ considered 
here, our definition practically coincides with the more general one given in pi]. 
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Ft 



ujjbd = jg- > 500 , (12) 

where F^q is the derivative of F w.r.t. (j> calculated at the present time. The latter bound does not depend on the 
particular cosmological trajectory at hand, since in contrast to (|ll|), it does not involve time derivatives, but only the 
present value of the field. 

Let us define now the form of F(tfi) that we consider here. We analyse non-minimally coupled scalar field models, 
defined as in |l4|. The term multiplying the curvature scalar R is made of a dominant contribution, which is a 
constant, plus one depending on (\>: 

^ )E i + ^ ) "^ o) ' (13) 

where 

F{ft>) = itf . (14) 

As it is evident from the previous equation, the present value of F(4>) has been setted in such a way that the 
gravitational constant reduces to the observed one; it is not our aim here to discuss how this result can be dynamically 
achieved. More general scalar-tensor theories in which Einstein gravity is an attractor have been discussed in the 
frame of Quintessence scenarios (see e.g. |27|). 

The second constraint ([l2]) turns out to be the dominant one for our models; as is very easy to verify, it gives rise 
to the following condition: 

10~ 2 

£<^7^- (15) 
2vG<po 

The only remaining quantity to specify in our model is the potential for the field </>, responsible for the the dominant 
part of the vacuum energy today. We adopt here inverse power-law potentials, as originally suggested by Ratra and 
Peebles R], and found in some phenomenological supersymmetry breaking models |f)|]3lf|, although cosine || and 
exponential || potentials have also been proposed in the literature. We take 

M A+a 

vto = — • ( 16 ) 

where the value of a > will be specified later and the mass-scale M is fixed by the level of energy contribution today 
from the Quintessence. 

Let us conclude this subsection by defining the quantities describing the Q energy density and pressure. By 
considering the kinetic and potential terms, we generalize in a straightforward way the definitions in ordinary General 
Relativity: 

P* = ^ + V{<t>) . (17) 
Similarly, the Quintessence pressure assumes the form 



4? - V{4>) . (is) 



However, as it is evident from equation (Q), there are other terms in the right-hand side that could be included into 
the Q energy density and pressure, coming directly from the coupling with R. Thus we define here the generalized Q 
energy density and pressure, that include all the NMC terms: 



±^ + V(cf>)-^f. ,1!). 



P* = iS$ > - v W + -2 + 0: r- (20) 

2cr a z a z 

We stress that, as we'll see below, in our model the NMC terms are negligible at the present time, so that, from the 
point of view of observational tests at low redshifts, there is no appreciable difference between the two definitions of 
energy density and pressure. 

In the next Section we will see how p<p,p<t, may differ substantially from p^^p^, while in general they become nearly 
identical in the MDE. 
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B. Perturbations 



As in 1 14 10], we adopt here the treatment of the scalar perturbations based on the formalism developed in |2q ] 
to describe the evolution of perturbations in the synchronous gauge. The perturbing tensor in Eq.(|^) is Fourier 
transformed and the generic component at wavevector k may be written as 



hij(k,r) = kikj/i(k,r) + (kikj - ^%)67?(k,r) , 



(21) 



where h denotes the trace of hij and r\ represents the traceless component; we omit the arguments (k, t) in the 
following. 

The perturbed Einstein equations read 



k 2 rj - -Hh = 
k 2 f] = 



a 2 Sp 
2~ ' 
a 2 (p + p)Q 



h + 2Hh - 2k 2 rj = -3a 2 Sp , 
h + 6ij + 2H(h + Grj) - 2k 2 r] = -3a 2 (p + p)cr . 

The perturbed density, pressure, velocity and shear in the synchronous gauge assume the form 

'p + 3p k 2 



Sp 



1 
F 



1 



HSF 



Sp f i uld + ^-f + -(-F 4> R + 2V (j> )S^-3 



SF 



Fh 



(22) 

(23) 

(24) 
(25) 



(26) 



Sp 



1 

F 



Spfluid 



\(F 4> R-2V 4 , 



a 2 ' 2--" 



SF HSF 



P- P 



2P 
3a 2 



SF 



1 Fh 
9^2" 



(27) 



/ \ i\ (P fluid + Pfluid)8 'fluid k 2 I — 

(p + p)o = j - 



-SF + HSF 



F 



(28) 



ip + P) a 



(P fluid + P fluid) O fluid 



2k 2 



F 



F ' 3^F\ SF + 3 ^ )+ 6 ) 



(29) 



where everything labeled with fluid contains contributions from all the species but Quintessence, and has the form 
as in ordinary General Relativity, obeying the perturbed continuity equations; we refer to JlO[ for a detailed definition 
of these terms. There remains the perturbed Klein-Gordon equation: 



5<f) + 2HS(f> 



k 2 + a 2 



-F A R + 2V 6 



6* = f + P,SR 



(30) 



Note the presence of the Ricci curvature scalar R in the left-hand side, as well as its perturbation SR in the right-hand 
side: the latter term is not trivial, since it contains S(f> (as it can be easily seen by using eqs. (|2^^) and (|3l]) below; 
from pq] we see that in the synchronous gauge it is in fact 



1 



SR = — [h - mh + 2k z r] 
so that the Klein Gordon equation, using Eq. (p4|) , assumes the complicated form below: 



(31) 



2H-^ 
2F 



2F lh 



HF d 



5<f)+ 



5 



-F^R + 



|| y-a 2 F^R - a 2 V^ + F^ 2 + F^ij> + HF^cf, ) + 



+F^ 



~6 ff V 



a 2 P - a 2 P , 2 2 
2 + 3* 



dpfiuid ~ ^ 



0ft . 



(32) 



This set of differential equations can be integrated once initial conditions on the metric and fluid perturbations 
are given; in this work we adopt adiabatic initial conditions for the various components as well as Gaussian and 
scale- invariant initial perturbations (see [ p~4|JTo) ] ) . 



III. TRACKING EXTENDED QUINTESSENCE 



Let us now turn to the study of tracking solutions in our model. As we already mentioned in the Introduction, the 
existence of this kind of trajectories was first pointed out by Ratra and Peebles 0. In Q , it was realized that they 
could be used in Quintessence scenarios to solve the fine-tuning problem on the initial conditions; this is in favor of 
Quintessence compared with the ordinary cosmological constant models, which remains affected by this problem. 

For what concerns the background parameters, we require the present closure density of Quintessence to be fi^ = 0.7, 
with Cold Dark Matter at VLcdm — 0.253, three families of massless neutrinos, baryon content = 0.047 and Hubble 
constant Hq = 100/iioo Km/sec/Mpc (or hioo — 0.65 thus keeping O(,/if 00 = 0.020). In section V we will explore 
the dependence of our results on the potential index a. Here we adopt a — 2. The case a — 1 was studied in p4[ , 
where however, we did not explore the properties of tracking solutions. As we will discuss in Section VI, however, the 
dynamical relevance of tracking solutions, becomes more and more relevant, as the potential slope increases, so, our 
main results of [jl4| are not essentially modified. 

In Figure fil we plot the evolution of the energy density of the cosmic fluid components (radiation, matter and 
Quintessence) as a function of redshift, in several models that we are going to describe now. The plotted curves 
are p r , p m and p§. The high, dotted lines represent radiation and matter, which do possess the well known scalings 
p r ~ 1/a 4 , p m ~ 1/a 3 . Equivalence occurs roughly at 1 + z eq ~ 5 x 10 3 . The other curves represent tracker solutions 
for the Quintessence energy density p^. We will describe in detail each curve in the next subsections; let us give here 
some general considerations. For each curve, the constant M in Eq.(p^) has been chosen to produce the required 
today. The NMC dimensionless coupling constant is taken as 

|£| = 1.5xl0- 2 , (33) 

so as to satisfy the experimental constraints dl2]), because in all the cases shown in Fig|l] the present value of <j) is 
4>o ~ 0.35Mp, where Mp = \j\fG. Also we study both positive and negative values of £. (f>o is reached starting from 

initial conditions 4>i, e g <C 4>o and ifibeg which can vary by several tens of orders of magnitude thanks to the capability 
of tracker solutions to remove any fine tuning of the initial conditions. In fact, it can be immediately noted from 
the curves in Fig.[y that, despite the huge range of initial values of the energy density, the Q component is going to 
dominate today at the chosen level. The final aspect we wish to point out before going to a detailed analysis of the 
dynamics in MDE and RDE, is that we plotted p$ in the figure, which is always positive, as it is evident from its 
definition (|l7]). This is not true for p^ defined in (|n]), because the NMC terms may be negative, and in fact they 
overcome the kinetic and potential energy, as we will see below; this feature was first found and discussed in []l4|| . 



A. TEQ-trajectories in the Radiation Dominated Era: i?-boost 



To understand what happens, let us study the Klein-Gordon equation during the Radiation Dominated Era (RDE). 
The first feature to note is the behavior of the Ricci scalar R. The dominant fluid is radiation, so that a ~ r, and 
this could yield the wrong idea that R = from Eq.(0). Actually this is not true as it is evident from Eq.(||): the 
first two terms Pfiuid — fluid take zero contribution irom radiation, but there is a residual contribution p m , from 
the subdominant matter component. Thus, there is a divergence R cx 1/a 3 as a — > 0, that holds until at least one 
particle species remains non-relativistic. Also, we will see in a moment that the other terms in Eq.(|3|), as well as the 
little dynamics implied by the overall factor 1/F do not change this argument; in conclusion the behavior of R in the 
RDE is 



G 



This implies that on the RHS of the Klein Gordon equation (||) the term multiplying R diverges as 1 /a and has the 
same sign of £ since <j) is assumed to be always positive. This generates a "gravitational" effective potential causing 
an enhancement of the dynamics of at early times, that we name "i?-boost" . The field accelerates until the friction 
term 2H4> reaches a value that is comparable to the i?-boost term on the RHS. After that, Q enters a phase of slow 
roll driven by the friction and the i?-boost terms only: 

2H^^fA, (35) 

The slow roll holds until the true potential energy from V becomes important in the Klein Gordon equation. This 
dynamics is manifest in Fig.||, where the absolute values of the four terms in the Klein Gordon equation are plotted. 
The thin solid line is the potential term, that is subdominant until 1 + z < 1000. The heavy dashed line is the friction 
term and the solid heavy one is the i?-boost term. In a very short time the friction grows from zero (the initial 
Q velocity is here taken to be zero) until it reaches the i?-boost term, setting the onset of the slow rolling regime. 
The thin dashed line is the Q acceleration, 0, which is positive and decreasing initially; then, it becomes negative 
(deceleration) at the cusp corresponding to 1 + z ~ 10 7 , and again positive at 1 + z ~ 10 3 when the potential energy 
becomes important and the tracker behavior in the MDE starts. 

It is quite simple to write an accurate analytic form for the solution during the i?-boost. By using Eqs.([54],[35|), and 

the behavior of a in the RDE, a ~ y/8irGp®./3 ■ t, we easily get 

<P = (t>be 9 exp [C(t - T beg )} , (36) 

where beg stands for the time when initial conditions are given, and 

c _ 3H$(tt C DM + Mb) £ ^ 
2y/8nGf#/3 

With our choice of parameters, C ~ 7 • 10~ 3 £ Mpc -1 , which implies that the exponent in ( |36| ) remains much smaller 
than one at all relevant times. Therefore, 

<f)~ fo eg [l + C{t - T beg )] . (38) 

The behavior expressed by Eqs.(^, |38|) corresponds to the dotted dashed line in FigJI], which mimics very tightly the 
i?-boost solution. Indeed in this phase the Q kinetic energy density dominates and Eq.(|3^) implies 

A 2 a z 2 a z 

The dot-dashed line in Fig.|l| has been obtained by inserting the value of C and (f>beg at the onset of the i?-boost. 

Note that the scaling as 1/a 2 is sensibly different from a truly kinetic dominated phase in minimally coupled models, 
for which we should have had 1/a 6 , from the continuity equation (Q) with — p<p- The reason why we have a different 
scaling of p$ is that the NMC terms in p$ are not negligible during all the RDE and the first part of the MDE, as we 
will show in a moment. Let us conclude the description of the i?-boost by estimating the time of its end. The latter 
occurs at a redshift z pot when the true potential energy becomes comparable with the kinetic energy (|39|): 

^ 2 beg C 2 (l + z pot ) 2 =V&(z pot )] . (40) 

Since the i?-boost is a very high redshift process, it covers a very tiny time interval, and in practice <f> does not move 
substantially from its initial condition (j)b eg during this phase. Therefore we can take V[4>{z pot )] ~ V{4>beg) in Eq.(|40"|). 
Moreover, for our inverse power-law potential we have V{4>b eg ) — V {4>o){4>Q 1 4>beg) 2 , and from the Friedmann equation 
today 87rGV(</>o)/3 = CI^Hq. Putting these ingredients together we can write an approximate formula for 1 + z pot 
giving the end of the i?-boost: 



Zpot ~ pnGC^tg ~ ' ( } 
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where the precise number has been obtained by making use of the actual value of the parameters in our case: it 
correctly corresponds to the i?-boost end in Fig.|l|. Note that the end of the i?-boost actually occurs when the 
Universe has already become matter dominated. 

Before concluding this subsection, let us return to the importance of the NMC terms in p^ in equation (p9[); as we 
already anticipated, we show now that they are dominant with respect to kinetic and potential energy densities until 
the first part of the MDE. Indeed, the NMC term in is —ZTLF^/a 2 and scales roughly as 1/a 3 , since Tt goes like 
1/a; on the contrary, the true potential V is roughly constant and the i?-boost kinetic energy scales like 1/a 2 ; as we 
explained. Because of these scalings, we expect that there exists a time tjssmc such that for r < Tnmc the NMC term 
is larger than both the kinetic and potential energy, while it becomes subdominant after this time. In fact, assuming 
that the time dependence of the field is as in the i?-boost solution (^), a simple calculation shows that the NMC 

energy density term —SHF^^/a 2 dominates p$ up to 



T NMC 



beg 



Q^PtoHq 



4>beg 



-i 1/3 



or 1 + z N mc — 850 



(42) 



that is long after matter radiation equivalence, as it is evident in Figj^. Note however that this is only an approxi- 
mation, since Eq.(o) tells us that the i?-boost solution is no longer satisfied at these redshifts; anyway this clearly 
shows that the NMT7 terms become subdominant only after the onset of the matter domination, so that the distinction 
between p^ and p<j, can be relaxed after tnmc- Our analysis is therefore approaching the MDE behavior of <fi, subject 
of the next subsection. 



B. TEQ-trajectories in the Matter Dominated Era 



Tracking solutions have been recently obtained in the context of Induced Gravity models pa] , whose only difference 
compared with NMC is the absence of a constant term multiplying the Ricci scalar in the gravitational sector of the 
Lagrangian, i.e. F(<j>) — £<f) 2 . We make two claims here. The first is that the same tracker solutions also exist in NMC 
models, provided in the MDE the scale factor can be expressed as a power of r. The second is that for F((j>) oc (fi 13 
the same solutions exist only if 8 = 2, as assumed so far. 

In fact, by looking at Eqs.(^, |), it is immediately realized that the difference between the IG and NMC models 
resides only in the 1/3F term multiplying the RHS in (0), since all the other terms involving F are derivatives with 
respect to either time or cj). Tracking solutions in |2q] have been obtained by assuming that the scale factor in the 
MDE is the square of the conformal time r, 



a(r) = a* — 



(43) 



(t* is a generic time) so that Eq.(|h_gets almost identically satisfied and disappears from the treatment. Indeed we 
show that in NMC models both Eq.M) and Ea.(M3) hold true, simply because matter is dominating and F{4>), playing 
the role of the gravitational constant m Eqs.(Q,|3}^ is slowly moving. Since 4>b eg <C 4>o — 0.35Mp, it is easy to see that 
the variation of the value of F throughout the whole cosmological trajectory is ^ 



Fbe 



Fo 



4.6 x 10" 



(44) 



To see this in another way, FigJ| shows the behavior of 7i, indicated on the vertical axis as (1 + zjH^ 1 , where 
H = at/ a, in three important moments of the cosmological evolution, namely present (top), decoupling (middle), and 
equivalence (bottom); the heavy, solid line represents TEQ with positive £, the thin solid line is TEQ with negative £, 
the short dashed line, instead, is ordinary Q, and the long dashed one is the cosmological constant, which is sensibly 
different from Q and TEQ since it is not dynamical. As it is evident, opposite signs of £ imply opposite behaviors for 
Ti, especially at small redshifts, compared with ordinary Q. At a given redshift, the shift in H. is due to the behavior 
of F((f>), which is less or more than 1/8ttG for positive or negative £, respectively, by an amount given by (E3). 



2 It is also useful to point out here that, because of the small change in the value of F(<f>) throughout the entire cosmological 
evolution, the well-known constraints on the variation of H at the epoch of nucleosynthesis are always satisfied in our models 
(see e.g. |l4| and references therein). 
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Therefore, Eq.([43|) holds with good accuracy, and all the solutions obtained in |2f| for IG models, as well as their 
stability properties, hold true also in our NMC case. 

Let us come now to our second claim. Let us assume for a while a form = ^ andF(0) = 1/(8ttG)+|(^-^), 
respectively. We are searching for power law solutions for the Quintessence 



(45) 



from (^) and (f45|), it is matter of simple algebra to check that Eq.(^) takes the form 



i o \ a 2 R ~, a , 9 M 4+a 

h 2 + 37) J = —m"- 1 + « 2 «^rr ■ ( 4 6) 



Also, from (|9|) the Ricci scalar in the MDE becomes R — 12/ (a • r) . It is easy to see that this is only satisfied if 

13 = 2 and a = - - 2 . (47) 

7 



a > and the field is positive, the following ac 



This proves our second claim. If conditions (J47) are satisfied, Eq.([4q) gives 0* in terms of 7,71, a; since in our models 



ditional condition is derived: 

7 2 + 3 7 - 12| > . (48) 
In this regime, the Quintessence energy density scales as 

P4> = P4>* — , e = — -z , (49) 
ya* / a + 2 

and its pressure is 

2 

= ~P<t> — ■ ( 50 ) 

For what concerns the stability of these solutions, we report here only a remarkable result (25|; solutions satisfying 
( [47| ) are stable under time dependent perturbations if and only if the following condition holds: 

-l-^<0, (51) 
a + 2 

hol ding for instance for any a > 0. In summary, for the class of solutions we are interested in, namely those satisfying 
( fi^,[46| ), the cases of interest here have been treated in |2f|; we have shown here that such solutions apply both to 
NMC and IG cases. 

Let us describe now in detail the solutions we see in Fig.jj]. Let us focus on the minimal coupling case first. The 
lowest, dashed line represents a minimally coupled case in which the behavior of the Q component during RDE is 
equivalent to a pure cosmological constant. The kinetic energy density is initially zero and remains largely subdominant 
with respect to the potential one. Thus cj> "freezes" for almost all the RDE, and leaves this condition only when its 
energy density reaches a fraction of about 1CP 3 of the critical one. Thereafter, Q joins the tracker solution in the 
MDE regime corresponding to 

7 = e = 1.5, (52) 

that satisfies the constraint ( f48| ) with f = 0. The substantial motion of the field from the initial condition to its final 
value occurs in this last phase. 

The thin long-dashed line represents a case in which the initial kinetic energy density is dominant with respect to 
the potential one, by 23 orders of magnitude. As it is evident, the field starts from an energy density comparable 
with the matter one. In this case the kinetic energy is redshifted away during the phase corresponding to the rapidly 
decreasing part of the trajectory in the figure; the scaling is easily found from the continuity equation, with ~ p^. 
Then the field freezes again before joining the tracker solution at 1 + z ~ 100. Note that because of this early stage 
of kinetic energy dominance, the field freezes at a value slightly greater than the initial condition 0t> eg = 10 -2 ; that 
is the reason why the flat part of the curve lies slightly below that corresponding to the previous case. Although not 
plotted in the figure, a roughly equivalent trajectory might have been obtained by requiring that the initial potential 
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energy density was comparable to the matter one. Note the very large set of initial energy values from which the field 
reaches the present state. 

Let us come now to the analysis of the tracker solutions in the NMC case. The heavy solid line has been obtained for 
£ = 1.5 x 10~ 2 . The RDE is dominated by the i?-boost: the field accelerates until the gravitational effective potential 
is reached by the cosmological friction term in the Klein-Gordon equation. The slow rolling sets in, and holds until 
the true potential becomes important, when the field freezes. After this early phase, p^ and p^, after having been 
much different in magnitude and sign, become indistinguishable, and join the tracker solution in the MDE reaching 
the required value today. 

The heavy dashed line corresponds to a case in which the initial kinetic energy density is larger by 23 orders of 
magnitude with respect to the potential one, thus starting from an energy amount comparable with the matter one. 
As for the Q case, in this condition the field undergoes an era of kinetic energy dominance until the latter is redshifted 
below the effective NMC potential energy and the i?-boost is set also in this case. The evolution of the field from this 
time on is the same as in the zero initial kinetic energy case. 

The dot-dashed line is in fact the i?-boost approximate solution (^6|). Note that it describes very well the i?-boost, 
and shows as the terms of order higher than the first in Eq.(|38|) become important only now, when the i?-boost is no 
longer active. 

The last point we want to stress, is that the heavy solid line is actually superposed with a solid thin line, which 
describes a case identical, but for negative coupling constant, £ = —1.5 x 10~ 2 . The reason why these two trajectories 
agree so tightly is that during the _R-boost the kinetic energy density, which dominates p^, is the same regardless of 
the sign of £, see Eqs. ( p6[ - p9[ ). The only slight difference between the two trajectories is when the i?-boost ends and 
the true potential starts to drive the dynamics of </>. For positive £, 4> is positive during the i?-boost and later, while, 
for negative £, </> is negative during the i?-boost. Although the curves describing p^ for positive and negative values 
of £ look very similar, this is not true for the perturbations they generate, as will be discussed in the next section. 



IV. COSMIC MICROWAVE BACKGROUND SPECTRA 



In this section we analyse the CMB temperature and polarization spectra of our TEQ model. 



In Fig. | 

perature 
from 



we plot the CMB spectra for the models discussed in the previous section; the top panel describes tem- 
uctuations; the bottom panel shows polarization spectra. The CMB spectrum coefficients are calculated 



C! =-!;: / ^|A„(fc,r )| 2 



f"' = 4;r / ^|A P ,(fc,T )| 5 



(53) 



where the quantities A^(fc, tq) and Api(k,To) are functions of the photon and baryon perturbed quantities, see e.g. 
ref. jl(J for detailed definitions. The spectra have been normalized at the COBE measurements at I = 10. 

The heavy, solid line in FigJ^ describes CMB spectra for the model corresponding to the same line in FigM The 
thin, solid line describes the same model, but for £ = —1.5 x 10~ 2 , again as the same line in Fig.pl The short-dashed 
line, in the left panels, represents a case of ordinary Quintessence with the same potential and 0^ = 0.7. The long- 
dashed line, in the right panels, describes a cosmological constant, A model, with J7a = 0.7. Before entering into 
the description of the various effects we find, let us stress that the effects of TEQ are quite large, with respect to 
both ordinary Q and A models. Also, the models plotted respect the experimental constraints discussed in Section II, 
having lujbd — 500 and Gt/G ~ 10~ 12 yr -1 . We want to mention that the effects are here considerably larger than 
what we obtained in our previous paper 14 1 ; the reason is twofold, in [ 14j we did not follow the tracker solution for 
the Q field, but for the shallower potentialwe considered there (a — 1) fnc size of the effects we find here would have 
been smaller anyway (see the discussion in Section VI). 

First of all, let us consider the dynamics of Q at low redshifts. By looking at Eq.(0), it is easily seen that at present 
Q obeys a sort of effective potential, caused by the true potential and by the curvauire-coupling term: 



V e ff((f> ) R) = V(<f>) 



F(4,)R 



(54) 



At low redshifts when Q starts to dominate, the dynamics of the Hubble length is suppressed since the universe is 
approaching a de Sitter phase. Therefore, by neglecting Ti in Eq.(^), we have R ~ 6H 2 . Also, from the Friedmann 
equation we get (87rG/3)V ~ fl^H 2 , and for inverse power-law potentials of the form dl6| ) we have dV/dcj) — —aV/4>. 
Thus the derivative of the effective potential (p4) takes the following approximate form: 



dV t 



eff 



dcj> 



{ 3o^0 



6£0 \H 



(55) 
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for positive £, both terms push toward increasing values of (j>, and we can immediately understand that the NMC 
term is comparable to the ordinary Quintessence one for 

IS " ' < 56 > 

Therefore, for £ = 1.5 X 1CP 2 and a — 2, we expect to have a 10% extra force coming from the NMC terms in Eq.fEt). 
Just the opposite happens if £ is negative, since the terms in ( f55| ) have opposite sign. Indeed what we see in Fig.gis 
that for opposite signs of £, the change in the low redshift dynamics of Hr 1 is also opposite. Also, in FigJ| it can be 
seen that at present roughly an order of magnitude separates the gravitational effective potential from the true one. 

The Integrated Sachs Wolfe effect (ISW, see e.g. B^] for a review) makes the CMB coefficients on large scales, 
or small £'s, change with the variation of the gravitational potential along the CMB photon trajectories. Since the 
gravitational potential is affected by the low redshift dynamics, we expect an increase or a decrease of the ISW for 
positive and negative £, respectively. Indeed this is precisely what we see in FigQ looking at the curves for £ < 10: 
positive or negative £ TEQ make the Ce larger or smaller than ordinary Q, respectively. 

The following simple calculation gives a good estimate of the amount of ISW in TEQ models. Take a cosmological 
scale comparable with the Hubble horizon today, so to be unaffected by acoustic oscillations. As it is known (see 
e.g. |32| and references therein), the ISW is essentially due to the change, between decoupling and now, of the 
gauge-invariant expression of the gravitational potential, \F In NMC theories ^dec is slightly different from the one 
calculated in minimally coupled theories, because it receives a contribution from the time variation of the gravitational 
constant, since ^ is proportional to G oc 1/F; going from decoupling to now, an increasing F (positive £) makes ^ 
decreasing, while the opposite is true for decreasing F. More precisely, in the limit \Fd ec — Fq\ <C Fq, \& changes by 
an amount 

The more the gravitational potential decreases, the stronger is the ISW power at £ < 10, which is roughly given by 
{ST/T)isw — 2(5\I/ [ p2| . Therefore, recalling that on superhorizon scales ST/T ~ on the CMB power spectrum 
the NMC contribution to the ISW can be estimated as 

6C t<io C ™io " °?<™ (ST/T)j smNMC - (6T/T)j sw F - F dec 2 

c i<w C^ 10 [ol / 1 ) ISW b Q 

the last equality has been obtained in our particular model, where at decoupling (j> <§; (f>o, so that 1— Fd ec /Fo ~ 8irG£(f)Q. 
Note again that, depending on the sign of £, the net effect can be an increase or a decrease of the CMB power. Indeed 
in Fig.0 we see that the above estimate is fairly well respected: since in the present case <po ~ 0.35Mp on large angular 
scales Eq.(j58|) predicts an effect of the order of 10%, as in the figure. 

Let us come now to the evaluation of the effects on the acoustic peaks in the CMB spectrum. There are essentially 
two effects, peaks amplitude and position, that we describe now. 

The amplitude changes in the opposite way for opposite sign of £. This feature can be understood as a normalization 
effect: we have seen that positive and negative £ give rise to increased and decreased ISW in TEQ with respect to 
Q models; correspondingly, when the spectra are normalized at £ = 10, a decrease and increase of the acoustic 
peak amplitude occurs. As a consequence, the magnitude of this effect is roughly at the level given by Eq. (|58|) . 
For completeness, we report also another mechanism that changes the acoustic peaks amplitude, which occurs at 
decoupling, instead of at low redshifts. As is evident in Fig.||, the size of the Hubble horizon at a given redshift is 
different in TEQ with respect to Q models. As we noted in p3, this implies that the horizon reentry for a given 
comoving scale is delayed for positive £ and anticipated for negative £; therefore, with respect to the ordinary tracking 
Quintessence, this implies a slight excess or deficit in CDM density when such scale is in horizon crossing, corresponding 
to a deficit and an excess in the radiation energy density. This slightly suppresses the acoustic oscillations in the first 
case, and enhances them in the second. 

Let us come now to the explanation of the acoustic peaks shift. The angular scales at which acoustic oscillations 
occur are directly proportional to the size of the CMB sound horizon at decoupling, that in comoving coordinates 
is roughly Tdec/v^i an d inversely proportional to the comoving distance covered by CMB photons from last scat- 
tering until observation, that is To — Td ec [p2[] . Therefore, since the spectrum multipoles scale as the inverse of the 
corresponding angular scale, we have 

£peaks OC — — . (59) 

Tdec 

In our models, tq and T dec shift essentially because of two reasons, namely the change in To due to the domination 
of the Q component today, and the behavior the Hubble length H^ 1 in the past. The first feature mainly makes 
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the Q and TEQ models different from a cosmological constant, since the latter does not have the kinetic degree of 
freedom. The second feature is simply related to the fact that, for a fixed value of the Hubble length today, in the 
past it possesses a behavior which is characteristic of the particular model at hand, because of the time evolution of 
the effective gravitational constant 1/F(4>) (see Fig. 3). Let us give a simple analytical formula for the acoustic peaks 
shift in TEQ models with respect to ordinary tracking Q. First, it can easily be seen that the conformal time r can 
be conveniently written as a function of the scale factor as follows: 

ad ~. (60) 
aa 

However, from the Friedmann equation (0), a scales as the inverse of \/F ; therefore, small changes SF <C F induce a 
change Six/ a ~ —(1/2)6 F/F, and consequently in the conformal time, which shifts by an amount 

sr=r r * , - Ta ir^, (6i) 

Jo aa[l - {1/2)SF/F} 2 J aa F V ' 



where we have defined the time change of F as follows: 

5F _ F{<f>) 



F 



(62) 



Of course SF/F at any given time depends on various details, but let us make the simplifying assumption that it is 
constant from the beginning until z — 2 and zero afterwards, empirically following what we see in Fig. 3. It is then 
immediate to deduce that Tdec changes as 



Instead, for what concerns tq we have 



1 5F 

Sr ^r z=2 --—, (64) 



because, according to our simplifying assumption, 5F/F w for z < 2. In conclusion, minding Eq.(|59|), and after 
some algebra, we get the shift in the acoustic peaks as a result of the time variation of the effective gravitational 
constant: 



A7 i te Q _ iQ , c p 

O^peaks peaks peaks 0" 

' Speaks " 2 F V To 



(65) 



Numerically we find t z= 2/to ~ 75%; also, we already mentioned that the change of the value of F during all the 
cosmological evolution can be written as 5F/F ~ — 87r£G0g. Therefore, for our specific model Eq.(|65|) becomes 

~ n£Gcf> 2 . (66) 

"peaks 

Note that, for our values of £ = ±1.5 x 10~ 2 , the above shift is at the level of ±6 x 10 -3 , which is in quite good 
agreement with the results plotted in the left panels of Fig. 4, that is ±5 x 1CP 3 . 

This completes our description of the TEQ features on the CMB angular power spectrum. We turn now to the 
analysis of what happens in the matter power spectrum today. 



V. MATTER POWER-SPECTRUM 



In Fig.g we plot the matter power spectrum for the same cases shown in Fig.[|. We can note immediately differences 
regarding both amplitude and turnover position. 

The A model has the highest spectrum. The main reason is the different growth of density perturbations || . In both 
Q and A models the perturbation growth is suppressed at low redshifts due to the domination of the vacuum energy, 
that tends to keep H constant therefore enhancing the cosmological friction in the perturbation equations, where 
almost everywhere the terms involving the first time derivatives of the perturbations appear multiplied by the Hubble 
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parameter. In Q and TEQ models this effect is considerably enhanced due to the magnitude of H which is greater 
than in A models at all redshifts, and in particular at the lowest ones, as it is evident from Fig.p|. Another independent 
cause that contributes to push the Q spectrum down with respect to the A model is the COBETnormalization. In fact, 
we have seen that the ISW effect is enhanced in Q models with respect to A ones, and the normalization to COBE 
implies subtraction of power to the true amplitude of the primordial cosmological perturbations. 

Let us come now to the difference between Q and TEQ models. The dynamics of the field at low redshifts is almost 
the same as in Q models, as it is evident again by looking at Fig|l]. Thus, the reason of the difference is to be searched 
for in the COBE normalization. Indeed, by looking at the low wavenumber region, which is the zone of non-processed 
scales, we see that the amplitude shift is roughly at 10% and —18% level for TEQ with respect to Q for negative and 
positive £, respectively; these numbers roughly agree with the ISW corrections that we estimated to come from the 
NMC terms in the previous section. 

Let us come now to the evaluation of the peaks shifts. As it is known (see for example the scale of the 
matter power spectrum turnover is essentially given by the scale entering the Hubble horizon at t : 
equivalence. The latter age is the same for all our models: 



le matter-radiation 



1 + Zeq = ^ a 5500 . (67) 

PrO 

However, we must take care of what was the Hubble horizon at the equivalence, since the Hubble radius follows 
different dynamics in the three cases that we are treating. In other words, the shift in the power spectrum turnover 
is given by 



Shurn ( SH 



" ... 



(68) 



In Fig.|| the different values of the Hubble horizon are displayed at the equivalence (bottom) . The Hubble horizon 
shift in Eq.(p8|) between the Q and A model is at the 18% level, which corresponds well to the power spectrum shift 
that we see m Fig.p, where Sktum/kturn — —20%. The same reasoning applies to explain the slight shift of the TEQ 
spectra turnover with resgect to the Q one. The Hubble horizon shift in Fig]| for positive £ is 2.6%, in good agreement 
with what we get in Fig.|| 2.3%; for negative £ we get 1.7% from both figures. 

We can give a rough analytical estimate of this effect by reasoning as follows. At the equivalence the Q ener 
density, both in Q and TEQ models, is negligible with respect to matter and radiation as it is evident in Fig.[y. Eq 
takes the form 

0,2 ^ P fluid (( . Q s 

H * 3F(f) ■ ( } 

Since in our model F(cf>) is smaller or larger than F(<fio) if £ is positive or negative, respectively, this implies that the 
Hubble length H^ 1 at the equivalence was different in TEQ models with respect to ordinary Q; the amount and the 
sign of the difference can be estimated roughly as follows: from Eq.(B^) we have 



5H 2 \ (28H\ 



= 1 - 8wG ■ F eq ((f> eq ) . (70) 



eq \ / eq 

Therefore, the matter power spectrum turnover shifts due to the NMC terms by the following amount: 

'6k turn \ _ fSH-^ 1 - 8TrGF eq ((t>) 



Hum J nmc V H / eq 



AttG&o , (71) 



simply by taking <f> eq <C 4>o as it is today. In our models, the quantity in Eq.(71) is roughly 2%, that agrees quite well 
with the numbers we obtain. 

For what concerns the shape of the power spectrum, there is no significant differences between Q and TEQ models, 
since this could arise only when the scalar field becomes important at low redshifts, where however the dynamics 
in the two scenarios is very similar. Notice that the oscillations that can be seen on the decreasing branch at high 
wavenumbers in figure || are just a residual of the acoustic oscillations of the photon-baryon fluid at decoupling. 

VI. VARIATIONS IN THE POTENTIAL SLOPE 

Varying the power a in the potential ( |l6|) implies varying the dynamics of according to Eq. ( fl9| ) . For increasing 
slopes of the potential, we expect that thelow redshift effects are enhanced correspondingly. 
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In FigJq we show tracking solutions for different exponents: dot-dashed line for a = 1, solid for a = 2, long-dashed 
for a. = 3. It can be immediately seen how the low redshift tracker branch possesses different slopes, according to 
Eq. (p9[) . Moreover, the i?-boost is the same for all the slopes, again according to the arguments made in Section II. 
But, ior increasing a, the i?-boost is abandoned earlier. The reason is the mechanism that stops the i?-boost itself, 
i.e. the fact that the Q potential energy starts to dominate. But, this happens earlier for larger exponents in Eq. (|l6|) , 
since the potential is steeper for (j> < <j)Q. 

Another important aspect that we must address is the Q equation of state, uu = p^/p^, since this is a parameter 
that is quite well constrained by the observations to be less than —0.6 today In fact, from Eq.(|50|) we see 

that small a make the equation of state similar to the pure cosmological constant case. The actual value that we 
find numerically is slightly different from what Eq.(|5^) would predict, because today the tracker regime has been 
abandoned, since Q has started to dominate the cosmic evolution. For a = 1,2,3, we find w<p = —0.75, —0.60 and 
—0.49, respectively, for positive £, and = —0.78, —0.66 —0.59, respectively, for negative £. 

Next, let us see what happens in the CMB. Fig.^ shows the behavior of the CMB temperature (top) and polarization 
(bottom) spectra for the same values of a as in Fig|| the left and right panels are for positive and negative £, 
respectively. Solid lines, both heavy and thin, refer to the same cases as in FigJ|; the dotted and dot-dashed lines 
describe models with a = 3 and 1, respectively. As we mentioned above, the low redshift ISW effect is enhanced by 
increasing a, merely because the low redshift dynamics is enhanced, as it is evident from Eq.(pq). This is a known 
effect since it does not involve directly NMC terms, but instead is a characteristic feature oTaynamical vacuum 
energy models flpj . The acoustic peaks region is suppressed correspondingly, by an amount that is roughly increasing 
linearly with a, due to the COBE normalization at low Ps. However, as it can be easily noted in FigM the strength 
of this effect is larger for positive £ compared with negative values. The reason is that negative £ values tend to 
decrease the ISW effect, as it is evident both in the effective potential ( p5| ) and in our estimate of Eq.(^8|). For what 
concerns the shift in the position of the acoustic peaks, the different potential mainly influences thelow redshift 
tracking behavior, thus affecting mostly To in ([39]); since large as enhance the low redshift dynamics, the peaks shift 
increases correspondingly, again as it is evident m Fig .J?]. 

In conclusion, the increased dynamics obtained by increasing a has the indirect effect of enhancing all the NMC 
effects. To see this, let us stress the part of the effects that come from genuine NMC terms. To this aim we build 
tracking solutions corresponding to those that we show in FigJq, but for the ordinary Q case, and calculate the CMB 



to quantify in detail the pure NMC effects in the spectra. These quantities, for the different values of a treated in this 
Section, are displayed as a function of the multipole I in FigM temperature (top) and polarization (bottom). The 
meaning of the curves is the same as in FigM The effects are larger for increasing values of a. In the upper box we 
can note the large ISW effect due to the NMC component only: it reaches roughly 30% by itself. Correspondingly, 
due to the normalization at small Ps, the acoustic peaks region is suppressed or enhanced for positive or negative £, 
respectively, even up to the 50% level. The enhancement of these effects, in particular the ISW one, for increasing 
a in in agreement with our formula (J5q), since for increasing a the field reaches larger values today: namely, for 
a = 1,2, 3, we have <fio = 0.2, 0.35 and TT5 Planck masses for positive £, and similar values for negative £. 

The same phenomenology obtained for the temperature anisotropies is found for the polarization (bottom panels); 
we will not plot the low Ps region, since even if the NMC effect is large, the absolute value of the polarization multipole 
moments is significantly high only around degree angular scales. 



In this section we will make a self consistent summary of the results obtained in this paper, focusing on the 
observable features of the cosmological models here investigated. 

The present paper is devoted to the study of theoretical aspects and observational imprints of cosmological scenarios 
in which a dynamical vacuum energy component dominates the cosmic evolution today, as suggested by high-redshift 
type la Supernovae data [Q. 

Extended Quintessence means the class of models for which the vacuum energy today is given by a scalar field 
which has an explicit coupling with the Ricci scalar R; we proposed this scenario in [ [14] and performed a preliminary 
study of its cosmological properties. In the present work we assume a Non-MinimaTCoupling (NMC) term of the 
form 



where £ is a dimensionless coupling constant, <fi is the Quintessence field, also indicated by Q, and 4>o is its value 
today. 




(72) 



VII. CONCLUSIONS 




(73) 
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As a first step, we investigated the cosmological trajectories in this scenario, focusing on the redshift evolution of 
the Q energy density. The Q potential is modeled as an inverse power law, namely V((j>) = M A+a / </>" , where a > 
and M is the energy scale chosen to have the amount of Q energy today. 

In the radiation dominated era we find that the divergence of R at early times generates an effective gravitational 
potential in the Klein Gordon equation, which drives the Q evolution in a sort of gravitational slow rolling. We named 
this era i?-boost. The dominant component in p^ during this epoch is the kinetic component. The i?-boost ends when 
this kinetic energy becomes comparable with the true potential energy: we have estimated that this time corresponds 
to a redshift 1 + z pot ~ 3000, i.e. the i?-boost lasts also for the earliest part of the matter dominated era. 

In the matter dominated era we obtained scaling solutions in NMC theories, which extend the validity of those 
studied in 123], for which p^ and the Q pressure p^, scale as follows: 



P<t> = Pqo 



\ a ) a + 2 a + 2 



This regime holds for almost all low redshifts, 1 + z < 100. At recent times, however, 1 + z < 3, the Quintessence 
component starts to dominate and the expansion accelerates. 

These results show that the Extended Quintessence is a tracker field, in the sense of the original works on 
Quintessence tracker solutions [0JM] ■ While the Quintessence potential energy density has to be chosen to produce 
the observed today, just as if Q was ordinary cosmological constant, the initial amount of p§ can vary over a very 
wide range, covering more than 20 orders of magnitude. 

We also investigate perturbation spectra in TEQ. The tracker dynamics imprints considerably larger effects com- 
pared with ordinary Tracking Quintessence and the cosmological constant, both on CMB and LSS. We gave quite 
accurate analytical formulas to describe these effects, that we summarize here. 

We have shown that in TEQ models the low redshift dynamics of the field is either increased or decreased depending 
on whether £ is positive or negative, respectively. This leaves its imprints on the large scale CMB angular power 
spectrum, where the signal is unprocessed by acoustic oscillations, and is affected by the dynamics of the gravitational 
potential (Integrated Sachs Wolfe effect, ISW). We derived the following simple formula describing this effect. If Ci 
is the coefficient describing the power of CMB anisotropics on large angular scales, I < 10, the shift in TEQ models 
compared with ordinary Q is given by 

<5Q<io 

12 (1 - 87rGF dec ) , (75) 



C, 



e<io 



where F is the NMC function in Eq. (|7|) . because of the tracking regime we generally have cj) beg -C cj) , so that the 
above equation assumes the simple form 



6C e<i 



, (76) 
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which means an effect going from 10 to 50% for the potential slope a increasing from 1 to 3. Note that for negative 
£ the ISW is inverted, a feature that is quite unusual in Quintessence models. 

The amplitude of the acoustic peaks varies because of the same ISW effect, if the overall power is normalized to 
COBE at £ = 10. There exists also another mechanism that changes the acoustic peaks amplitude, which occurs at 
decoupling rather than at low redshifts. The Hubble horizon value at a given redshift is different in Extended vs. 
ordinary Quintessence models, as already noticed in fll4j . This is because in order to reach the same value today 
the Hubble parameter followed a different dynamics inthe past; we have shown here that the horizon reentry for a 
given comoving scale is delayed for positive £ and anticipated for negative £. This implies a slight suppression of the 
acoustic oscillations in the first case, and an enhancement in the second, as a consequence of the different amounts of 
CDM when the baryon and photon densities are oscillating. 

TEQ imprints also affect the position of the acoustic peaks. The angular multipoles at which CMB acoustic 
oscillations are occurring are inversely proportional to the comoving size of the CMB sound horizon at decoupling, 
proportional to Trf ec , and directly proportional to the comoving distance between last scattering and us, tq — Td ec - 
Both these quantities arc affected by the TEQ dynamics. The CMB sound horizon is proportional to the Hubble 
horizon at decoupling, which, as we just mentioned is different in the various models. The comoving distance between 
last scattering and us is also affected by the TEQ low-redshift dynamics. We derived the following simple analytical 
formula giving the shift in the positions of the acoustic peaks, which accurately describes our numerical results: 

61 ^GF{4> dec ) - 1 2 

j g = ^G^o • (77) 

Note that again, depending on the sign of £, Eq.(Q) may imply a peak shift towards smaller or larger angular scales. 
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For what concerns the LSS effects today, the first point to make is that since at the COBE scales, where we normalize 
our models, the CMB spectrum is affected by the ISW, the amplitude of perturbations on LSS scales is suppressed if 
£ > and increased if £ < 0, by an amount that is approximately the opposite of the quantity in Eq.(|76|). Second, the 
turnover scale position k tU rn is shifted, being inversely proportional to the Hubble horizon at the equivalence. The 
size of this effect can be simply estimated by the Friedmann equation: 



Again, this is roughly AttG^q because, thanks to the tracking behavior, <f) eq -C <fio- 

This completes the list of our main results. A very final comment we want to make is that CMB anisotropics 
have been shown here to be considerably affected by the dynamical properties of the vacuum energy and by the Non- 
Minimal Coupling terms present in the Lagrangian. This is an important novelty in the literature on Quintessence 
models, which deserves as much attention as the dependence of the CMB spectra on the kind of dark matter, value of 
the primordial perturbation spectral index, etc.. This is made possible by the fact that the CMB itself is more than 
a snapshot of the early universe; its properties in fact are also determined by a sort of line-of-sight integration over 
a very long part of the cosmic evolution, which is then able to tell us about the time variation of the gravitational 
constant and the dynamics of the vacuum energy. 

With higher and higher accuracy achieved in CMB data, going from the encouraging:, present-day ones [^8|,|9| to the 
formidable performance of the Microwave Anisotropy Probe and Planck mission |17[], we will be able to gerexciting 
new insights into the structure of gravity as well as the nature of the vacuum energy. 




NMC 



1 - 8irGF c 
2 



(78) 
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FIG. 1. Matter, radiation, Quintessence energy densities as functions of redshift. Dotted lines: matter and radiation; 
heavy solid: Tracking Extended Quintessence; thin dashed: tracking minimally coupled Quintessence; heavy dashed: initially 
highly energetic Tracking Extended Quintessence; thin long-dashed initially highly energetic tracking Quintessence; dot-dashed: 
i?-boost analytic. 
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FIG. 2. Absolute value of the Klein Gordon equation terms: heavy solid: gravitational effective potential; heavy dashed: 
friction; solid: potential; dashed: acceleration. 
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FIG. 3. Snapshots of conformal Hubble horizon evolution, heavy solid: Tracking Extended Quintessence with positive £; 
short dashed: tracking minimally coupled Quintessence; thin solid: Tracking Extended Quintessence with negative £; long 
dashed: cosmological constant. Top: present; middle: decoupling; bottom: equivalence. 
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FIG. 4. CMB angular power spectra. Top: temperature; bottom polarization. Heavy solid: Tracking Extended Quintessence, 
positive £; thin solid: Tracking Extended Quintessence, negative £; short dashed: minimally coupled tracking Quintessence; 
long dashed: cosmological constant. 
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FIG. 5. Present matter power spectra. Heavy solid: Tracking Extended Quintessence, positive £; thin solid: Tracking 
Extended Quintessence, negative £; short dashed: minimally coupled tracking Quintessence; long dashed: cosmological constant. 
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FIG. 6. Matter, radiation and Quintessence energy densities for different potential slopes. Long dashed: matter; short dashed: 
radiation; heavy solid: Tracking Extended Quintessence, a = 2; heavy dotted dashed: Tracking Extended Quintessence, a = 1; 
heavy long dashed: Tracking Extended Quintessence, a = 3. 
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FIG. 7. CMB angular power spectra for different potential slopes. Left (right): positive (negative) £. Solid: Tracking 
Extended Quintessence, a = 2; dotted dashed: Tracking Extended Quintessence, a — 1; short dashed: Tracking Extended 
Quintessence, a — 3. 
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FIG. 8. CMB angular power spectra: relative difference between Tracking Extended Quintessence and ordinary tracking 
Quintessence spectra (SCi/Ci = Cj E ® /Cf — 1) as a function of the inverse power of <j> in the potential. Top (bottom): 
temperature (polarization); thick (thin) lines for positive (negative) £. Solid: extended vs. ordinary tracking Quintessence, 
a = 2; dotted dashed for a = 1, short dashed for a = 3. 
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